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Overview 


• Flight Dynamics Mechanisms -Ballistic and Powered Flight 

• Ballistic Flight with Nonlinear Dynamics (“Chaos”) 

- 2-Body, 3-Body, 4-Body 

• Powered Flight 

- High Thrust (Apollo) vs. Low Thrust (SMART-1) 

• Creative Combinations 

• Trajectory simulations 

- Low thrust trajectory from GTO to Moon 

- Gravity Assist 

- Weak Stability Boundary 

• Mission Trade-offs 

• Conclusion and Future Work 
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Ballistic 2- and 3-Body Motion 


• 2-Body (Keplerian) motion: 

- Classical ellipses 

- Ellipse evolves if central body is non-spherical like Earth or lumpy like the Moon 

• 3-Body Motion: Earth-Moon-S/C 

- Gravity assist 

• Switch from Earth-S/C to Moon-S/C and back to change orbit energy and 
angular momentum (esp. inclination) 

• Double Lunar Swingby acrobatics allows us to explore translunar space, and 
rotate orbit axis to track Sun (far left image) 

- Unstable Motion around Collinear Libration Points LI and L2 

• Unstable dynamics allows small change to have large long-term effects 

• Low energy transfers from Earth to Moon via EM LI “corridor” (Conley, 1968, 
center image) 

• First LPO mission ISEE-3 (1978) designed by Farquhar et al. at NASA GSFC 

• Stable and unstable manifolds can be used for low-energy transfer, e.g. by 
ARTEMIS mission (right image) 
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Ballistic 2- and 3-Body Motion (cont’d) 

• Earth-Moon Resonance 

- Highly Elliptical Earth orbit resonant with Moon to get regular small changes from 
Moon 

• Stable: Keep apogee far from Moon to gain long-term stability and eliminate 
maintenance maneuvers (IBEX, Left image) 

• Unstable: put apogee near Moon and let Moon raise apogee (e.g. SMART-1 
and trajectory concept in Right Image) 

- This type of asymptotic trajectory was described by Deprit & Henrard (1 965) 



IBEX, launched from a Pegasus 



Figure 6.12: An Earth- to-Moon transfer departing from a LEO. This solution 
requires a cost equal to Av = 3244 m/s and a time of transfer At = 65.1 days. 
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Ballistic 4-Body Motion 


• 4-Body Motion: Sun-Earth-Moon-S/C 


- Weak Stability Boundary (WSB) transfers 

• Switch from Earth-Moon-S/C to Sun-Earth-S/C and back to change orbit 
angular momentum (orbit shape and inclination) 

• S/C is sent to vicinity of Sun-Earth LI or L2, roughly 1 .5 million km from 


Earth 


• Developed by Belbruno and Miller at JPL 

• First demonstrated with Japanese Hiten mission (1990) 

- Primary mission in Highly Eccentric Orbit with apogee 476,000 km to 
swing past Moon. 

- Belbruno proposed lunar swingby to send S/C near Sun-Earth LI , then 
return to Moon with no delta-V required for lunar capture. (Left image) 

• In 2011 , WSB was used to transfer GRAIL S/C to lunar orbit (Right image) 



NuneftcaJ illation of a bnllkUic capiare transfer trajectory 
for the Japanese spacecraft Hiten: ediplic plwr projection, stiffs 
direction indicated at Earth injection. < Ikttwu Miller | IWfi 
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ARTEMIS Transfer Trajectory - PI 


• Use of Multi-body Environment to raise periapsis to lunar orbit 

• Two lunar gravity assist separated by 13 days for flip of apoapsis direction 

• Deterministic Deep Space Maneuver (DSM1) was performed 33 days later 

• All maneuvers target the Earth-Moon libration insertion state 

• Unstable Lissajous manifold towards the Earth-Moon system 

• An l_ 2 Lissajous insertion orbit maneuver of 2.56 m/s 



‘Back-flip’ 


Lunar Flyby #1 
31 -Jan-2010 
@ 12600 km altitude 


Rotating Sun-Earth Coordinates: 
jn-plane view 


Maximum Range 
(6-Jun-20l0). 

1 50 million km 

To Sun 

◄ 


Trajectory Correction 
Maneuvers (TCMs) 

• 0.31 m/s on 2 -Feb-2010 

• 0.18 m/s on 20-Apr-2010 

• 0.18 m/s on 20-Jun-2010 

• 0.65 m/s on 19-Jul-2010 

• 2.24 m/s on 18-Aug-2010 


Deep Space Maneuver (DSM) 
• 7.31 m/s on 10-Mar-2010 
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Total Time from ORM-1 to LL2 

Insertion 387 Days 

Total Time from 1 st Fly-by to LL2 

Insertion: 204 Days 

ELI and EL2 are about 1 .5 

million km from Earth 


lunar Flyby M2 
13-1 


Out-of-Plane 


LL2 Lissajous 
Orbit Insertion 
23-Aug-20l0 
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ARTEMIS PI TCM Unstable Manifolds 


• Unstable manifold and the effect of TCM5 on the return orbit 

• Original planned and the corrected post TCM5 trajectory shown 

• For the actual maneuver, the difference is slight in terms overall design but the shift in the 
general direction of the flow is consistent with a new manifold 

• The post TCM5 path with a small shift in direction guaranteed that PI would reach its Cartesian 
goal at the proper epoch 
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Targeting Chaos using Resonance 


Bollt & Meiss (1995) proposed a low-energy 
transfer to the Moon using chaotic behavior 
near Earth-Moon LI point. 

• Transfer would take two years. 

Schroer & Ott (1997) proposed low-energy 
trajectories exploiting resonances that reduced 
the transfer time to 377 days and 293 days. 

• Center image shows Highly Eccentric Orbit near 
5:3 then 3:2 resonance with Moon 

Bernelli-Zazzera et al. (2003) found a simpler 
trajectory that used resonance to reach the 
Moon in 61 days (Bottom image). 

NOTE: These trajectories are in the Planar 
Circular Restricted 3-Body Problem. Three- 
dimensional motion and the Moon’s orbit 
eccentricity present additional challenges 
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FIG. 6 . Controlled chaotic trajectories to a stable lunar orbit, (a) Starting 
from a circular orbit 59 669 km above the earth's center (time required for 
transfer: 377.5 days), (b) Starting from the period three orbit denoted by 
open circles in Figure 5 {time required for transfer: 293 days). 



f igure 6.12; Ail KarLfa-lo-Moon Srwitfw departing from * I.Ett Thi» Hulnlfc.n 
reigiirt?* h m<4 rep ml Us At 1 = 3244 m/a mill a I slur i:J IrsuisfiT if = Ki.l dttlta. 
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Powered Flight 


• High Thrust Rockets and Engines, Apollo style 

- Uses Hohmann Transfer to minimize impulsive 
velocity change (delta-V) 

- Requires delta-V to escape Earth orbit, and to 
enter Moon orbit 

- Allows fast transfer to the Moon, within 3 to 4 days 

• Low Thrust Engines, SMART-1 style 

• Continuous thrusting for long intervals, spiraling 
out from Earth/in toward Moon 

• Typically takes months 

• Does not require carrying large propellant mass, 
so can deliver more mass to final destination 



APOLLO 11 FLIGHT PROFILE 



* For best results, we combine ballistic flight 
with powered flight options 
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Combining Low Thrust, Resonance and 

Stable Manifolds 


• Belbruno’s 1987 paper on Lunar 
Capture Orbits combined low 
thrust with use of LI corridor. (Top 
image) 

• Mingotti, Toputto, Bernelli-Zazzera 
(2006) found a trajectory that starts 
with low thrust from GTO, reaches 
a resonance orbit, then rides the 
corridor around Earth-Moon LI to 
an L2 orbit (Bottom image) 

• Martin and Conway (2010) 
employed a similar approach to 
reach lunar orbit 
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SMART-1 


Fig. 3. From long eclipse to moon capture (Thrust: dotted line; Diamond: resonance; Apogee: plus sign). 

Trajectory design "based originally on interior transfer" 
described in Belbruno's 1987 Lunar Capture Orbits paper 

Launched 27 Sept 2003 into GTO from French Guiana (low 
inclination) 

Continuous thrust to raise perigee out of radiation belts, using 
SEP 

Thrust around perigee to raise apogee 

Use Moon resonances to achieve lunar capture (5:1, 4:1, 3:1, 

2 : 1 ) 

Spiral into moon orbit 

Arrived in lunar orbit Feb 2005 (13 month transfer) 

' ' -500000-400000-300000-200000-100000 0 100000 200000 300000 400000 500000 

X coordinate (km) (direction of moon oscending node w.r.t. EQ2000) 




Fig. 4. From long eclipse to moon capture (Thrust: dotted line; Diamond: resonance; Apogee: plus sign). 
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SMART-1 -like Approach 


After examining the various options, we decided to focus on a trajectory similar to 
SMART-1 

We have a software tool "LTOC" developed by JAQAR Engineering, that specifically 
models low-thrust transfer from GTO to lunar orbit 

The LTOC solution has the following segment: 

- Raise apogee using low-thrust arcs near perigee 

- Spiral out from Earth using continuous low-thrust along velocity vector 

- Coast to time encounter with Moon 

- Lunar capture 

- Spiral in toward Moon using continuous low-thrust along anti-velocity vector 

In the following charts we present some results from that LTOC solution 

We also used the LTOC solution as a starting point to create a trajectory in 
STK/Astrogator, which uses an gravity assist to achieve plane change for lunar 
encounter 

NOTE: These are trajectory concepts, not yet fully developed trajectory design. 
Propulsion system design and constraint analysis (e.g. power, communications) 
remain to be done. 
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LTOC Trajectory 



Earth centered edge view 
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Sunlight Conditions 



U**gM 'tfTCGi 


LTOC design does not 
consider Sun lighting 
constraints, but it is an 
important consideration for 
Solar Electric Propulsion 
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STK Scenario: WSB with Lunar Gravity Assist 
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STK Scenario: WSB Edge View 
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Trade-off Table 



Impulsive 

(mono-prop) 

76 

220/ 

750 

2.9 

30 % 

* Single impulsive to achieve a 
Weak Stability Boundary 
transfer 

Low Thrust 
Continuous - 

314 

2500 / 3500 

1.5 

15 % 

• Uses Lunar assist for orbit 
plane change 

GTO 






Low Thrust 
-GTO/ 
JAQAR Tool 

428 

2500 / 3880 

2 

20% 


Low Thrust 
periapsis 
maneuvers 

1048 

2500/1500 

0.6 

6 % 

• Minimal shadows 

• Uses Lunar assist for orbit 
plane change 


-GTO 
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Conclusions and Future Work 


• We have examined the diverse mechanisms of ballistic and powered flight to send 
a spacecraft from GTO to lunar orbit 

• We chose to focus on a SMART-1 like trajectory, but that does not exclude the 
possibility of using a different approach. 

• The LTOC tool by JAQAR Engineering allowed us to develop an initial trajectory 
from GTO to polar lunar orbit. However LTOC is limited in its design capability 

• We then developed a similar scenario in STK, using a lunar gravity assist to achieve 
plane change (like Hiten). Thus far we have trajectory concepts but not a trajectory 
solution 

• We have begun a study to examine trade-offs in transfer duration, mass to mission 
orbit, range from Earth, radiation exposure 

• For further trajectory design we need to define the propulsion system 

• We need to develop our software tools to accurately model the many spirals 
around Earth during low thrust orbit raising 

• We would also like to add dynamical systems analysis to the trajectory design, to 
identify useful manifolds 
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Backup Charts 
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Transfer Trajectory - PI 


• Use of Multi-body Environment to raise periapsis to lunar orbit 

• Two lunar gravity assist separated by 1 3 days for flip of apoapsis direction 

• Deterministic Deep Space Maneuver (DSM1) was performed 33 days later 

• All maneuvers target the Earth-Moon libration insertion state 

• Unstable Lissajous manifold towards the Earth-Moon system 

• An l _2 Lissajous insertion orbit maneuver of 2.56 m/s 



Back-flip’ 


Lunar Flyby #1 
31 -Jan-2010 
@ 12600 km altitude 


Maximum Range 
(6-Jun-20l0) 

1 50 million km 


To Sun 

◄ 


Deep Space Maneuver (DSM) 
• 7.31 m/s on 10-Mar-2010 

Lunar Cubes 2013 


Total Time from ORM-1 to LL2 

Insertion 387 Days 

Total Time from I s1 Fly-by to LL2 

Insertion: 204 Days 

ELI and EL2 are about 1 .5 

million km from Earth 


lunar Flyby #2 
13-1 


Out-of-Plane 


Rotating Sun-Earth Coordinates: 
In-olane view 


Trajectory Correction 
Maneuvers (TCMs) 

• 0.31 m/s on 2 -Feb-2010 

• 0.18 m/s on 20-Apr-2010 

• 0.18 m/s on 20-Jun-2010 

• 0.65 m/s on 19-Jul-2010 

• 2.24 m/s on 18-Aug-2010 


LL2 Lissajous 
Orbit Insertion 
23-Aug-20l0 
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Transfer Trajectory - P2 


Use of Multi-body Environment to raise periapsis to lunar orbit 
Single lunar gravity assist 

Three Deterministic Deep Space Maneuvers were performed 
All maneuvers target the Earth-Moon libration insertion state 
Unstable Lissajous manifold towards the Earth-Moon system 
An L 1 Lissajous “insertion” orbit maneuver of 0.25 m/s 


Rotating Sun-Earth Coordinates: 



Out-of-Plane 


In-plane view 


To Sun 

◄ 


LL1 Lissajous 
Orbit Insertion 
22-Oct-2010 


Lunar Fly-by 
28-Mar-2010 @ 8070 km 
altitude 


Total Time from ORM-1 to LL1 
Insertion 458 Days 
Total Time from Fly-by to LL1 
Insertion 208 Days 
ELI and EL2 are about 1 5 
million km from Earth 


Trajectory Correction 
Maneuvers (TCMs) 

• 2.15 m/s on 20-Jul-2010 

• 0.72 m/s on 2-Aug-2010 

• 0.31 m/s on 1 -Oct-2010 

• 0.25 m/s on 12-0ct-2010 


Deep-Space Maneuvers (DSMs) 

• 3.68 m/s on 13-May*2010 

• 23.3 m/s on 1-Jun-2010 

• 3.28 m/s on 9-Sep-2010 


Maximum Range (18- 
Jun-2010) 
t 2 million km 
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PI Trajectory Correction Maneuvers 


• As the transfer trajectory was flown, Trajectory Correction Maneuvers (TCM) 
were required to adjust for previous maneuver performance 

S Function of temperatures and general propulsion system performance 
S Spacecraft attitude errors and spin pulsing directions (pulse width) 

S Navigation errors 

• To find a maneuver location, a DC 
and an optimization was run to determine 
the optimal maneuver location 

S Determine an a priori maneuver 
location 

S Achieve an intuitive feel for 
maneuver results 
S Meet DSN coverage 
S Minimize DV 
S Optimization includes S/C maneuver direction constraints 

S Analysis reflects manifold behavior (up DV on one side can be down DV on the other) 

• Final planning was a combination of Optimization, Manifold Information, and 
experience in order to meet spacecraft constraints and insertion target goals 
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PI TCM Stable Manifolds 


• Consider only the outbound arc and TCM5 maneuver. 

• From the DSM, PI follows the original outbound path to the location of TCM5 

• TCM5 will shift the spacecraft to a different path, one that can be envisioned in terms of a different 
(orange) manifold 

• Subsequent to the DSM, and along the outbound trajectory two outbound manifold arcs emerge 
from the TCM5 location 

• These two manifolds represent the potential outcomes from (1) flow along the optimal path and (2) 
the alternative that incorporates a possible TCM5 maneuver 



PI Planned Stable Sun-Earth Lunar Cu !F¥ IWand Post TCM5 Stable Sun-Earth Manifold 23 
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PI TCM Unstable Manifolds 


• Unstable manifold and the effect of TCM5 on the return orbit 

• Original planned and the corrected post TCM5 trajectory shown 

• For the actual maneuver, the difference is slight in terms overall design but the shift in the 
general direction of the flow is consistent with a new manifold 

• The post TCM5 path with a small shift in direction guaranteed that PI would reach its Cartesian 
goal at the proper epoch 
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LunarCube Low-Thrust Trajectory Design 

TO Calculator Optimisation Output 


AFir s t Cut 


Earth-Moon Low-thrust Optimisation 

Solution does not comply or only partially complies with user inputs 


Using an optimization low thrust simulation tool 
Identified a sample design 


Launch from Earth Orbit - free orbit and free epoch 

Launch Mass: 12 kg 

Earth orbit at TLI: 

TLI time: 19 Sep 2016 13:19:58 

TLI delta-V: 2.72393 km/s 

semi-maj. axis: 27477.976 km 

Eccentricity: 0.76061 

Pericenter: 199.707 km 

Apocenter: 41999.972 km 

Inclination: 27.199 deg 

RAAN: 19.572 deg 


0 


s 


□ZD 


tTT^ctory in MEE2000 



349570 

329807 

310043 

290280 

270516 

250753 

230989 

211226 

191462 

171699 

151935 

132172 

112408 

92645 

72881 

53118 

33354 

13591 

-6173 

-25936 

-45700 

-65463 

-85227 

-104990 

-124754 

-144517 

-164281 

-184044 

-203808 

-223572 

-243335 

-263099 

-282862 

-302626 

-322389 

-342153 

-361916 

-381680 

-401443 


-396054 -334239 -272425 -210610 -148795 


-25166 36649 98463 160278 222093 283907 345722 407537 469351 531166 5929 : 


Arg. of peri: 72.909 deg 

True anomaly: -92.192 deg 

Transfer orbit: 

TLI time: 19 Sep 2016 13:19:58 

Phase 1: Apo raising (thrusting between TA -101.38 and +101.38 deg) 

Phase ends: 27 Jan 2017 2:21:8 

Phase 2: s.m.a. raising (continuous tangential thrusting) 

Phase ends: 22 Sep 2017 23: 7:16 

Phase 3: Targeting to Moon (thrusting between TA -101.38 and +101.38 deg) 
Phase ends: 24 Sep 2017 6:26:49 

Phase 4: Coast phase (no thrusting) 

Phase ends: 25 Sep 2017 6:26:49 (entry Moon SOI) 

Phase 5: LOI (tangential thrusting within Moon frame) 

Transfer time: 468.96047 days 

LOI time: 1 Jan 201812:23:2 

LOI delta-V: 1.16050 km/s 

Moon orbit parameters: 
semi-major axis: 2953.200 km 
Eccentricity: 0.37756 

Perilune: 100.000 km 

Apolune: 2330.000 km 

Inclination: 90.000 deg 

RAAN: 224.467 deg 

Arg. of peri: 270.000 deg 

Lat. of peri: -90.000 deg 

B-plane B*R: 56344.359 km 

B-plane B*T: 53016.530 km 
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Total delta-V: 3.88443 km/s 
Final Mass: 10.000 kg 
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LunarCube Low-Thrust Trajectory Design 






Simulation uses a 10 kg mass and a 1.4 mN thruster with an Isp of 2500 sec 
Thrusting is continuous along the velocity vector until it achieves a fairly circular 
orbit just below the lunar orbit 

The optimizer permits the s/c to enter into the a multi-body dynamical region to use 
the lunar gravity to 'rendezvous' with the moon at one of the relative node crossings 
Optimal approach LT maneuver back propagates to an optimal lunar distance to 
begin the low thrust lunar targeting sequence to capture into an elliptical orbit that 
is further decreased in semi-major axis by continuous thrusting after lunar capture 
The final orbit achieved is 


Keplerian - Moon Fixed 

sma: 9705.0226067659096000 km 

ecc: 0.8076956451338013 

inc: 109.2596992424557 deg 

C3: -0.5051818295180461 km A 2/sec A 2 

Rad. Peri: 1866.3181113559920000 km 

Rad. Apo: 17543.7271021758280000 km 

Periapsis is ~ 128 km. The period is 1-day 


RAAN: 104.5707969374672 deg 

w: 5.582162902317398 deg 

TA: 359.9999999991558 deg 

Energy: -0.2525909147590231 km A 2/sec A 2 

Vel. Peri: 2.1791742423983500 km/sec 

Vel. Apo: 0.231822595774652 km/sec 


The total transfer time of 314 days, from Sept 19, 2016 to July 30, 2017 

Total DV required is 3.458 km/s and the fuel mass required is 1.316 kg (13% of initial 

wet mass) 
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LunarCube Low-Thrust Trajectory Design 


• Earth Departure and Lunar Capture 
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LunarCube Low-Thrust Trajectory Design 



Altitude 


Satellite-Lunar Cubesat- 20 Mar 201 3 09:24:45 


4.0e5 — 


3.5eB -- 


3.0eB — 


2.5eB 


2.0eB 


1.5eS — 


I.OeB 


0.5eB -- 



May Jun Jul Aug 


Altitude (km) 
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LunarCube Low-Thrust Trajectory Design 


Shadows 


Umbra Eclipse Duration statistics 
Min Duration 3 Mar 2017 21:16:35.078 
Max Duration 4 Feb 2017 08:38:19.742 
Mean Duration 

Sateliite-Luna r_ Cub esa t: Lighting Times - 20 Mar 2013 09:27:55 


3 Mar 2017 21:34:11.342 

4 Feb 2017 09:22:47.661 


17.604 minutes 
44.465 minutes 
36.215 minutes 


Umbra 


Penumbra 


\ \ t 


Sunlight 


2016 


+ 

Oct Nov Dec Jan 2017 Feb Mar Apr 

Time (UTCG) 


May 


Jun 


Jul 


Aug 


Sunlight (UTCG) 


LunETiTiM T ^3 


Umbra (UTCG) 
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LunarCube Low-Thrust Trajectory Design 


• Simulation uses a 10 kg mass and a 1.4 mN thruster with an Isp of 2500 sec 

• Thrusting is performed only near periapsis along the velocity vector until it achieves 
a lunar encounter at apoapsis 

• The optimizer permits the s/c to enter into the a multi-body dynamical region to use 
the lunar gravity to 'rendezvous' with the moon at one of the relative node crossings 

• Optimal approach LT maneuver back propagates to an optimal lunar distance to 
begin the low thrust lunar targeting sequence to capture into an elliptical orbit. 

• Periapsis is ~ 3200 km. The period is 4-day. Orbit is weakly captured and will impact 
after several days. 

• The total transfer time of 1048 days, from Sept 19, 2016 to Aug 3, 2019 

• Total DV required is ~1.5 km/s and the fuel mass required is ~0.6 kg (6% of initial 
wet mass) 
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LunarCube Low-Thrust Trajectory Design 

• Sept 2016 to Jan 2019 


• l hrust only at periapsis^ from ta H'd to l30 cfeg 

• Three distinct oeriapsis maneuver scenarios 





LunarCube Low-Thrust Trajectory Design 
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LunarCube Low-Thrust Trajectory Design 


Satellite-Lunar_Cubesat _perigeeDV_Pianechange_Capture - 26 Mar 2013 11:31:50 



R_Mag {km} 



LunarCube Low-Thrust Trajectory Design 

Shadows 


Umbra Times 


Start Time (UTCG) Stop Time (UTCG) Duration (hr) 


19 Sep 2016 13:19:58.000 19 Sep 2016 13:19:58.000 0.000 

23 Aug 2019 10:22:13.855 23 Aug 2019 10:50:47.199 0.476 

28 Aug 2019 00:11:30.889 28 Aug 2019 00:35:30.555 0.400 

31 Aug 2019 12:07:23.194 31 Aug 2019 12:21:19.196 0.232 


Global Statistics 


Min Duration 19 Sep 2016 13:19:58.000 
Max Duration 23 Aug 2019 10:22:13.855 
Mean Duration 
Total Duration 


19 Sep 2016 13:19:58.000 
23 Aug 2019 10:50:47.199 


0.000 

0.476 

0.277 

1.108 
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LunarCube Low-Thrust Trajectory Design 

Earth Departure Dv at Periapsis ~ 720 m/s 


System 


Universe Mission Parameters Orbit Files Uiews <Last> 


76.2095 


days 



7b . 2895 4«yi 


WSB Transfer 


Sun-Earth Line 


Solar Rotating Coordinates 
Side view 


Earth 


Moon 


CubeSat 


Lunar orbit 


Solar Rotating Coordinates 
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